Microscopic two-nucleon overlaps and knockout reactions from 12 C 
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The nuclear structure dependence of direct reactions that remove a pair of like or unlike nucleons 
from a fast 12 C projectile beam are considered. Specifically, we study the differences in the two- 
nucleon correlations present and the predicted removal cross sections when using p-shell shell-model 
and multi-fttj no-core shell-model (NCSM) descriptions of the two-nucleon overlaps for the transitions 
to the mass A=10 projectile residues. The NCSM calculations use modern chiral two-nucleon and 
three-nucleon (NN+3N) interactions. The np-removal cross sections to low-lying T=0, 1() B final 
states are enhanced when using the NCSM two-nucleon amplitudes. The calculated absolute and 
relative partial cross sections to the low energy 10 B final states show a significant sensitivity to 
the interactions used, suggesting that assessments of the overlap functions for these transitions 
and confirmations of their structure could be made using final-state-exclusive measurements of the 
np-removal cross sections and the associated momentum distributions of the forward travelling 
projectile-like residues. 

PACS numbers: 24.50.+g, 25.70.Mn, 21.10.Pc 
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I. INTRODUCTION 

Properties of the wave functions of pairs of nucleons 
in a mass A + 2 projectile can be probed using sudden 
two-nucleon removal reactions, that exploit fast, surface- 
grazing collisions of the projectile with a light target nu- 
cleus. The sensitivity is to the wave functions of the 
nucleons at and near the surface of the projectile. When 
combined with 7-dccay spectroscopy, partial cross sec- 
tions of populated final-states in the mass A reaction 
residues can be determined. This direct reaction mecha- 
nism, its cross sections, and their distributions with the 
momenta of the forward travelling residues, are now be- 
ing exploited as a spectroscopic tool in studies of the evo- 
lution of nucleon single-particle structure near the Fermi 
surfaces of some of the most exotic nuclei; see for ex- 
ample Refs. [Ml and citations therein. The reaction 
observables used are, currently, inclusive with respect 
to the final states of both the removed nucleons and of 
the struck light target nucleus. More exclusive measure- 
ments, e.g. of the light charged fragments in the final- 
state, may in the future provide additional probes of the 
projectile structure [||. 

A detailed discussion of the two-nucleon removal reac- 
tion mechanism, its eikonal reaction-dynamical descrip- 
tion, and the cross sections and their momentum distri- 
butions, in the case of high-energy two-nucleon removal 
from 12 C, was presented in Ref. [||. A feature of this 
model description is that the removal cross sections in- 
volve only elastic interactions of the projectile residues 
with the target but sums of contributions from both elas- 
tic and inelastic interactions of one or both nucleons with 
the target @. New data, for the sd-shell nucleus 28 Mg 
and the 28 Mg(— 2p) reaction, have quantified these differ- 



ent contributions experimentally Q, and have confirmed 
that the relative importance of these different processes 
to the cross sections are consistent with the predictions 
of the eikonal dynamical model. This has provided an 
important additional test of the reaction model. In the 
earlier work for 12 C [f|, the theoretical comparisons used 
the sums of these removal contributions and p-shcll (Ohu 
shell-model) structure calculations were used to construct 
the required ( 10 X(JJ, T)| 12 C) two-nucleon overlaps. The 
PJT § and the WBP @ shell-model effective-interaction 
Hamiltonians were used. 

Key elements of that analysis are also relevant here: 

(i) The reaction is geometrically selective [l(| and the 
two-nucleon removal cross sections will be enhanced if 
the projectile ground-state has components with pairs of 
nucleons with strong spatial correlations (localization). 

(ii) The available experimental cross section data, from 
high-energy primary-beam measurements, are inclusive 
wit h respect to the populated bound states of the residues 
following np, nn and pp removal [ill [l2l | . The data, at 
three energies, reveal a significant enhancement of the 
ratio of unlike-pair yields, <7_ n p, to those for the like- 
nucleon pairs, o~— nn and C—pp. This enhancement is sig- 
nificantly greater than that expected based simply on the 
numbers of available 2N-pair combinations (i.e. a factor 
of 8/3). For example, the experimental <7- np : cr_ nn ratio 
was 8.54 for the data set with a 12 C beam of energy 2.1 
GeV per nucleon. Some (but not all) of this enhance- 
ment could be explained as due to the pair-correlations 
already generated in Ofouj p-shell-model overlap functions 
and because a larger fraction of the nn-removal strength 
leads to unbound 10 C final states. However, the experi- 
mental (J-np remained factors of 1.45 to 2.2 larger than 
the theoretical model calculations for the three available 
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data sets [TT|, E3 • Table HI reproduced from Ref. Q , 
shows both the p-shell-model results and data. 

(iii) The shapes and widths of the reaction residues' 
momentum distributions have both diagnostic and spec- 
troscopic value, being indicative of the total angular mo- 
mentum, /, the total orbital angular momentum L, and 
hence, with (LS)I coupling, also the total spin S carried 
by the removed nucleon pair fl4j . 

(iv) The calculated cross sections for the T—l states 
common to all three residues, namely the first + and 
2 + states, are essentially equal. Minor differences in the 
calculations arise from the small differences in the em- 
pirical separation energies for each system. Unlike the 
np- removal case, these calculated inclusive T=l cross 
sections were reasonably consistent with and were frac- 
tionally larger than the data values for o~- nn and cr_ pp , 
suggesting that the deficit in the theoretical cross sec- 
tions in the np channel reflects, primarily, a failure of 
our description of the overlaps ( 10 B(JJ,T = 0)| 12 C) for 

the transitions to the T=0, 10 B final states. 

Here we exploit the eikonal reaction model in the 
isospin formalism 0, [l(| HH for the removal of the like 
(T=l) and unlike (T=0,1) pairs of nucleons from 12 C. 
In particular, we will investigate the effect on reaction 
observablcs when using ab-initio multi-?kj no-core shell- 
model (NCSM) descriptions for the two-nucleon overlaps. 
We contrast these with the earlier Ohui p-shell shell-model 
results of Ref. [f|. We discuss the np and nn removal 
channels to 10 B and 10 C. These channels share the same 
T = 1, 0^ and 2+ final states, see e.g. Fig. [TJ whereas 
the pp channel contains two additional T — l states. 

In Section [II] we reiterate some specific features of 
the reaction description for 12 C projectiles. The nec- 
essary formalism has been presented elsewhere @, [l(| 
to which readers are referred. We follow the notations 
used in these earlier works. The chiral effective field the- 
ory (EFT) two- and three-nucleon (NN+3N) interactions 
and the microscopic no-core shell-model (NCSM) calcu- 
lations used to construct the improved overlap functions 
will be discussed in section lnTl These include calculations 
in which the chiral 3N interaction is switched off in an at- 
tempt to understand the impact on observables of these 
3N interaction terms in the starting Hamiltonian. The 
new results and predictions for 12 C reaction observables 
are analyzed in section IIVI and a summary is presented 
in section [V] 



II. CARBON-INDUCED REACTIONS 



tainty and were taken at high energies where the eikonal 
model description of the reaction dynamics is most reli- 
able. These data were obtained for reactions of a car- 
bon beam and carbon target at 250, 1050 and 2100 
MeV per nucleon incident energies. As mentioned above, 
these data (and related data for other light projectile nu- 
clei) show a significant enhancement in their T=0,1, np- 
rcmoval production cross sections over those with T=l, 
the nn and pp removal cases (see e.g. Table |T|. This 
observed enhancement is of particular interest as a po- 
tential signal and a measure of strong np-correlations at 
the nuclear surface. 

The primary motivation for the present study is the 
implementation and first assessment of improved micro- 
scopic descriptions of the two-nucleon overlap functions 
and their implications for rap-correlations and the calcu- 
lated reaction yields and observables. The relevant fi- 
nal states of the A=10 residual nuclei and their spins 
and isospins are shown in Fig. [1] The known spectrum 
of low-lying states in 10 B also contains several negative 
parity states. These arc not expected to be populated 
in the two-nucleon removal reaction mechanism. More 
exclusive measurements would be needed to confirm this 
expectation. 

The isospin formalism developed in Refs. [f|[l(| is used 
here. The description of the reaction and the parameters 
used for the 12 C and mass A=10 residue densities, etc. 
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Consideration of two-nucleon removal from 12 C is valu- 
able given the availability of both conventional shell- 
model and ab-initio NCSM structure descriptions. The 
residual nuclei 10 C, 10 B (and 10 Be) were also extensively 
studied and so establish a valuable benchmark. In ad- 
dition, the existing experimental cross section measure- 
ments [ll|, [llj], although inclusive with respect to the 
residue final states, have relatively small quoted uncer- 



FIG. 1. States of the mass A = 10 residues populated in 
the two-nucleon removal reactions considered. The spin and 
isospin labels (J/,T) are indicated. All states included are 
of positive parity. Levels assumed to be part of the isospin 
multiplet are connected by dashed lines. The lowest parti- 
cle thresholds are also indicated. States above the a-particle 
threshold in 1 B are expected to decay via a emission, with 
the exception of the 5.164 MeV, T=l, JJ = 2+ state which 
has an 84% 7-decay branch. 
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TABLE I. Calculated and experimental cross sections for two-nucleon removals from 1 C, for projectile energies of 250, 1050 
and 2100 MeV per nucleon, from 0. All cross sections are in mb. The calculations use p-shell shell-model wave functions from 
the wbp effective interaction. 



Energy 




10 Be 






10 c 






io B 




MeV/u 


C-2JV 


&exp 


Cexp/o'-2N 


C-2N 




~exp 1 l O -2N 


C-2N 


G exp 


&exp/o'-2N 


250 [11] 


7.48 


5.88±9.70 


0.79±1.30 


5.80 


5.33±0.81 


0.92±0.14 


21.57 


47.50±2.42 


2.20±0.11 


1050 [12] 


6.62 


5.30±0.30 


0.80±0.05 


5.13 


4.44±0.24 


0.87±0.05 


19.27 


27.90±2.20 


1.45±0.11 


2100 [12] 


6.52 


5.81±0.29 


0.89±0.04 


5.04 


4.11±0.22 


0.82±0.04 


19.02 


35.10±3.40 


1.84±0.18 



are the same as were discussed and tabulated in Ref. |5| . 
The key approximations were as follows: (a) The removal 
of nucleons is sudden, their co-ordinates frozen during 
the short timescale of the relativistic, surface-grazing re- 
actions, (b) The no-recoil (heavy residue) approxima- 
tion is made @ . The inclusion of core recoil affects only 
the diffractivc-stripping terms of the reaction mechanism 
whose contributions are a significantly smaller fraction of 
the cross section («30%) than the two-nucleon stripping 
terms, for which recoil is not an issue. Thus, recoil can 
play, at most, a minor role on computed cross sections. 
It is, in any case, the fractional changes of the cross sec- 
tions between the NCSM and the p-shell calculations, and 
between T=0 and T=l configurations that are of most 
interest and significance to this work. 

The primary difference here is a now considerably- 
extended set of two-nucleon amplitudes (TNA) C^ T aris- 
ing from the multi-shell set of available 2N-configurations 
a = [f3i, /y. Here the index j3 = (nlj) denotes the spher- 
ical quantum numbers of each active single-particle state 
in the model space. We evaluate the cross sections for 
transitions from the projectile initial (ground) state i, 
with spin (Jj,.Mj), to particular residue final states /, 
with (Jf,Mf). The all-important two-nucleon overlap 
function for removal of two nucleons 1 and 2 is 

- E CFViMJfMflJiMi) 

IfiTa 

{TTT f Tj\T lTl ) [# 1 (l)®V/s 2 (2)]F; • (1) 

Here Jf=0 + and thus Jf=I, the 2N total angular mo- 
mentum. We discuss later the use of Woods-Saxon or 
harmonic oscillator radial functions for the single-particle 
orbitals ipp. The 2N correlations under discussion arise 
in Eq. (fT]) from: (a) somewhat trivially, the antisym- 
metry and angular momentum coupling of nucleon pairs, 
and (b) the possible coherent pair-enhancement arising 
from the weights and phases of the TNA that contribute 
to each Jf final state. The latter and their sensitivity to: 
(i) the NN+3N effective interactions used, and (ii) the 
NCSM model-space dimensions, are studied here. 

Importantly, unlike the earlier p-shell-only analysis, 
the NCSM wave functions and overlaps now include ac- 
tive single-particle orbitals of both parities, including up 
to N max = 6 major oscillator shells. The angular correla- 
tion function of the two nucleons thus contains both even 



and odd Legendre polynomial terms in the angular sep- 
arations of the two nucleons [3, [l|| with the likelihood 
of a greater spatial localization of pairs at the nuclear 
surface. The degree of such increased correlations is dis- 
cussed further in the results section where we also show 
the two-nucleon joint-position probabilities in the im pac t 
parameter plane, Vj. (si, s 2 ), as were defined in Ref. [14j . 
These display the two-nucleon spatial correlations as de- 
livered by the projectiles and as are seen by the target 
nucleus that induces the reaction. 

Nucleon removal may occur via either elastic (diffrac- 
tion) or inelastic (stripping) interactions of nucleons with 
the target nucleus, the former leaving the target in its 
ground state. The latter lead to cross section contribu- 
tions that are inclusive with respect to all other target 
final states. Two-nucleon removal events can involve: (a) 
both nucleons making inelastic collisions with the tar- 
get, (b) there being one inelastic and one elastic col- 
lision, or (c) both nucleons suffering clastic collisions. 
Events (b), referred to as diffraction-stripping, are iden- 
tified in the reaction's absorption cross section but re- 
quire a projection-off bound states for the elastically in- 
teracting nucleon (for details see Refs. @,[nl). When 
truncated basis, single major shell, shell-model calcula- 
tions have been used to generate the TNA, all active 
singlc-particlc orbits in the shell were included in this 
projection operator. Here, when using the NCSM cal- 
culations in a very large basis, we have limited the pro- 
jection to the Op-shell orbitals, as being the appropriate 
bound states set. Purely elastic 2N- removal events (c) 
were estimated, as previously @, from the stripping and 
diffractivc-stripping cross sections. Typically, they con- 
tribute to the cross sections at the level of < 5%. 

The absorptive eikonal S-matrices, that largely deter- 
mine the volume of the two-nucleon overlap function that 
is sampled in the reaction, were calculated by folding the 
target, nucleon, and reaction residue point- nucleon den- 
sities with a zero-range effective nuclcon-nucleon inter- 
action. Further details may be found in Ref. Our 
primary interest is in the np removal cross sections. Since 
both the 0^ and 2^ T = 1 states are populated in both 
10 B and 10 C we may also extract the cross section for the 
nn channel. Very minor differences in these partial cross 
sections will arise from (a) the use of a 10 B rather than a 
10 C S-matrix, and (b) from the small binding energy dif- 
ferences. The latter are not accounted for in the present 
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calculations where harmonic oscillator radial wave func- 
tions are used. 



III. NN+3N INTERACTIONS AND OVERLAPS 

The p-shell (OHlj) shell-model calculations and overlaps 
were described in Ref. [f| and were computed using the 
code OXBASH For the present work, a series of no- 
core-shcll-modcl calculations, each for a given number of 
major oscillator shells, N max =0, 2, 4 or 6, were carried 
out using two chiral EFT NN+3N interaction choices, 
denoted NCSMl and NCSM2 in the following. 

The calculations used interactions derived within the 
chiral effective field theory (EFT) approach. In particu- 
lar, the chiral N 3 LO NN interaction of Ref. [HIHl was 
used with or without the chiral N 2 LO 3N interaction [2l| 
in the local form of Ref. [22|. These interactions were 
softened by the similarity renormalization group (SRG) 
technique [23j-|25| , where a unitary transformation is used 
to suppress the off-diagonal matrix elements (controlled 
by a parameter A). The SRG interaction induces higher- 
body interaction terms. These induced terms were kept 
up to the three-body level. It has been shown (2(|[27| that 
four- and higher-body terms are negligible for light nu- 
clei although some evidence for four-body induced terms 
was observed in 12 C calculations with one of the inter- 
actions used here (NCSM2) [13]. In NCSMl the NN+3N 
Hamiltonian used a 3N cutoff of 400 MeV and used pa- 
rameters fitted to the 3 H lifetime and the 4 He binding 
energy [H] . In NCSM2 the 3N cutoff was 500 MeV and 
the parameters were fitted to the lifetime and binding 
energy of 3 H [29| . In both cases the SRG was carried out 
using A=1.7 fin -1 , although the NCSM2 calculations were 
also performed with A=1.88 fm _1 to verify the SRG-A 
independence, i.e., to confirm the unitarity of the SRG 
transformation. The subsequent NCSM calculations used 
an harmonic oscillator (HO) basis with an angular fre- 
quency hio=16 MeV. The mass-dependent parameteriza- 
tions of the oscillator frequency Huj = 45^4 _1 / 3 — 25A -2 / 3 , 
agreeing with charge radius observations, suggest a value 
of fta>=14.9 MeV [31|, in reasonable agreement with the 
value used here. 

In the case of the NCSM2 parameterization, the calcu- 
lations were also repeated, and denoted as ncsm3, when 
the chiral 3N interaction in the starting Hamiltonian 
was switched off, but with the SRG-induced 3N effects 
(with A=1.7 fin -1 ) included. Again, the HO frequency 
of hio=l6 MeV was employed. Using these NCSM3 TNA 
we can make a first assessment of the impact on calcu- 
lations/observables of the inclusion, or not, of the chiral 
3N interaction in the starting Hamiltonian. 

It should be noted that the 10 B structure poses a par- 
ticular challenge to ab initio calculations. In particular, it 
had been observed that standard accurate NN potentials 
predict incorrectly the ground-state of 10 B to be l + 0, in- 
stead of the experimental 3 + 0. The present calculations 
with the chiral N 3 LO NN potential (NCSM3) suffer from 




FIG. 2. (Color online) Experimental excitation energies of 
1 B are compared to the different calculations used in the 
present work: chiral N 3 LO NN (ncsm3), chiral N 3 LO NN 
plus N 2 LO 3N with the cutoff of 400 MeV (ncsmI), and chiral 
N 3 LO NN plus N 2 LO 3N with the cutoff of 500 MeV (ncsm2). 
The -/V max =6 space was used in calculations shown in the first 
four columns. The SRG A parameter is indicated. The HO 
frequency of M7=16 MeV was used in all calculations. 



the same problem. Only after including the chiral N 2 LO 
NNN term, with the 3N cutoff of 500 MeV NCSM2, docs 
one get the correct ground state spin. Interestingly, the 
weaker chiral N 2 LO NNN with the 3N cutoff 400 MeV, 
NCSMl, fails to invert the l + and 3 + states, also pre- 
dicting the wrong 10 B ground state spin. Sec Fig. [2] for a 
comparison of 10 B excitation energies from different cal- 
culations used in this paper. Also in the figure, the sta- 
bility of the spectra with respect to the SRG A variation 
and the size of the model space V max is demonstrated 
for the NCSM2 case. The situation is somewhat reversed 
in 12 C, where the Hamiltonian NCSM2 with the stronger 
3N interaction over-binds 12 C by several MeV and over- 
corrects the splitting of the l + and 4 + states 27j. Us- 
ing the weaker 3N interaction (NCSMl) both the binding 
energy and excitation energy description improves. Fur- 
thermore, this Hamiltonian (NCSMl) also describes the 
binding energies of oxygen and calcium isotopes [28j very 
well. The stronger 3N interaction NCSM2, on the other 
hand, provides a very good description of lighter nuclei 
(A < 10), resolving even long-standing analysing power 
problems in p— 4 He scattering [3(|. These observations 
suggest that our knowledge of the 3N interaction in par- 
ticular is incomplete and additional terms, such as those 
at the N 3 LO of the chiral perturbation theory, must be 
included. Further, the mass region of ^4=10 — 12 is ideal 
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to test the details of the nuclear Hamiltonians. 



IV. RESULTS 

Reaction calculations were carried out using the ex- 
tended TNA sets derived from the two NN+3N effective 
interactions outlined in section Hill for each of N max =0, 
2, 4 and 6, and for the nn and np-removal channels. In 
Table [TT] we first show the calculated inclusive cross sec- 
tions for the p-shell wbp and (the most complete) NCSM 
N max =6 calculations at the three energies of the avail- 
able data. The NCSM calculations use harmonic oscilla- 
tor radial form factors with hui=16 MeV, as used in the 
NCSM basis. The experimental inclusive cross sections 
were shown in Table HI We note that when using the P JT 
p-shcll effective interaction the inclusive 10 B yield at 2.1 
GeV per nuclcon was 18.73 mb as compared to the 
19.02 mb shown in Table [TT] for the wbp case. Thus the 
two p-shell calculations are entirely consistent in this in- 
clusive cross section observable. We note that the p-shell- 
model calculations, taken from use Woods-Saxon ra- 
dial wave functions when constructing the two-nucleon 
overlap functions. If hu!=16 MeV oscillator functions are 
used the calculated cross section is reduced by ~ 7% for 
the lowest energy 3 + , and by < 3% for higher energy 
states. The smaller changes for the higher residue ex- 
citations are due to the absence of any binding energy 
effect when a single oscillator frequency is used. 

The cross sections, now exclusive with respect to the 
A=10 final states, are shown in Table Hill for the highest 
energy of 2.1 GeV per nucleon. The two NN+3N NCSM 
cases (with N max =6) and the WBP case are shown to- 
gether with the NCSM3 case (i.e. without the chiral 3N in- 
teraction in the starting Hamiltonian). In rap-removal to 
10 B the cross sections are shown for the six positive par- 
ity gamma-decaying final states below the first nucleon 
threshold. However, the first 2 + , T=l state is known 
to decay by a-cmission with an I a = 16% branch. This 
branching has been accounted for in the inclusive (T- 2 n 
values presented. The cross sections for population of 
the three higher lying T=0 states (see Fig. [1]) are not 
included since these states are reported to decay by a- 
emission (with I a = 100%). It was assumed that these 
states do not contribute to the 10 B yields. 

Table IIIII reveals significant sensitivity of the yields of 
the low-lying 10 B states to the interactions assumed. The 
ratio of the T=0, if to 3 + ground-state yields is reversed 
between NCSMl and NCSM2. The absence of the 3N in- 
teraction in the NCSM2 starting Hamiltonian, the NCSM3 
case, leads to a quite significantly enhanced ratio of the 
T=0 if to 3 + ground-state cross sections when com- 
pared to the full NCSM2 Hamiltonian. To confront these 
detailed model predictions requires more exclusive mea- 
surements with good statistics. 

The six contributing 10 B partial cross sections are also 
plotted in Figs. [3] El and [B] Fig. Elshows the calculations, 
using NCSMl, for N max =0, 2, 4 and 6, and the stabi- 



lization and the essential convergence of the calculated 
partial cross sections (upper panel) and the full-width at 
half-maximum widths of their momentum distributions 
(lower panel) with increasing N max . Based on this ob- 
served convergence, seen for all of the NCSM cases, we 
have presented, in the main, only the final N max =6 re- 
sults. Fig. @] shows N max = 6 NCSM2 calculations for 
two different SRC A values, A=1.7 fm -1 and A=1.88 
fm -1 . The essential independence of the cross sections 
and the FWHM on the SRG-A confirms the unitarity of 
the SRC transformation, i.e., we are really investigating 
predictions of the initial chiral interactions. Based on the 
results of Ref. [27|, the SRC A dependence is expected 
to be even weaker for the NCSMl and NCSM3 cases. Fig. 
[5] shows the calculations from the NCSMl, NCSM2 and 
NCSM3 Hamiltonians, all with N max =6. The significant 
variations predicted in the widths of the momentum dis- 
tributions for the different final states and the sensitivity 
of the T=0, 10 B final state yields to these effective inter- 
actions, are evident. The momentum distribution widths 
from the NCSM interactions are broadly consistent, and 
also with those of the conventional truncated-basis p- 
shell model interactions (shown in Fig. [5]). They are 
systematically wider than the latter by rj 25 McV/c that 
we attribute to the incorrect asymptotic of the oscillator 
wave functions. These similarities reflect the relative in- 
sensitivity of the different LS'-fractions in the overlap to 
the interactions used. The notable exception is the sec- 
ond T = 0, 1 + excited state, when the width using the 
NCSM interactions is somewhat wider than that from the 
conventional shell-model interactions. The comparison of 
the NN+3N and p-shcll shell-model results are shown in 
Fig. M 

The final-state inclusive cross sections for the NCSM 
interactions are all larger than those from the truncated- 
space WBP interaction. All three NCSM interactions re- 
sult in similar summed cross sections to the lowest two 
T = states, between 16 and 17 mb, but the distribution 
of this strength between the two states shows significant 
variations. The cross sections to the T = 1, + state is 
also larger when using the NCSM interactions with vari- 
ations between the predictions of the interactions, thus 
reinforcing the need for final-state-exclusive cross section 
measurements. The detailed nature of the sensitivity of 
the TNA, and hence the final-state branching ratios and 
cross sections, to the details of the interactions is com- 
plex, but precise measurements could provide a path to 
probe this sensitivity and toward constraining the under- 
lying interactions. 

The ratios of the theoretical model and experimental 
cross sections, R S (2N) = a exv la t h are R s (2n) = 0.66 and 
R s (np) = 1.40 for the NCSMl interaction and R s (2n) = 
0.56 and R s (np) = 1.35 for NCSM2. Calculations for 
two-nucleon removal in exotic sd-shell isotopes typically 
overestimate the experimental observations by a factor of 
two with R S (2N) = 0.5 [(J. As mentioned above, we do 
not include core recoil and ccntcr-of-mass effects for the 
structure amplitudes, so do not draw specific conclusions 
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TABLE II. Calculated inclusive cross sections for two nucleon removals from 12 C, for projectile energies of 250, 1050 and 2100 
MeV per nucleon. All cross sections are in mb. The TNAs used were calculated using the p-shell and wbp interaction, and the 
NSCMl and NSCM2 NN+3N interactions with N max =Q>. The experimental cross sections were tabulated in Table|T] 



Energy 




io c 






io B 




MeV/u 


WBP 
< J -2N 


NCSMl 
"-2N 


NCSM2 
"-2N 


WBP 
0"-2iV 


NCSMl 
°"-2JV 


NCSM2 
°"-2JV 


250 


5.80 


7.10 


8.48 


21.57 


28.19 


29.91 


1050 


5.13 


6.31 
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FIG. 3. (Color online) Calculated exclusive cross sections (up- 
per panel) and full width at half maximum (FWHM) widths 
of the momentum distributions (lower panel) for the ground 
and 7-decaying final states of the 10 B residues, following np- 
removal at 2100 MeV per nucleon. The TNAs used were ob- 
tained using the NCSMl NN+3N starting Hamiltonian (dashed 
lines). Calculations are for N max =0 (open triangles), 2 (red 
triangles), 4 (open circles) and 6 (blue circles). 



FIG. 4. (Color online) Calculated exclusive cross sections (up- 
per panel) and FWHM widths of the momentum distributions 
(lower panel) for the ground and 7-decaying final states of the 
B residues, following np-removal at 2100 MeV per nucleon. 
The TNAs used were obtained using the NCSM2 NN+3N start- 
ing Hamiltonian (dashed lines). Calculations are for N ma x=6 
with the SRG A=1.7 frn -1 (open circles) and A=1.88 fm _1 
(red circles). 



about these absolute values of R S (2N). However, both 
effects should be independent of the reaction channel, 
and the former will not affect the strongest mechanism, 
i.e. 2N-stripping. 

To consider the impact of the significantly larger 
model spaces introduced using the NCSM amplitudes, 
it is meaningful to consider the cross section ratios 
^■m* / a -2N M 1 i- e -> t ne relative enhancement of a par- 



ticular channel when moving from the truncated basis 
wbp to the NCSM interactions. For the NCSMl interac- 
tion, this ratio is 0.81 for the nn channel and 0.76 for the 
np channel. The corresponding numbers for the NCSM2 
interaction are 0.69 and 0.73. Evidently, the use of the 
larger basis amplitudes set enhances the cross section. 

A significant fraction of these differences results from 
changes in the p-shell amplitudes, as are shown in Ta- 
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FIG. 5. (Color online) Calculated exclusive cross sections 
(upper panel) and FWHM widths of the momentum distribu- 
tions (lower panel) for the ground and 7-decaying final states 
of the 10 B residues, following np-removal at 2100 MeV per nu- 
cleon. The TNAs used were obtained using the NCSMl (blue 
circles) and NCSM2 (open circles) NN+3N starting Hamilto- 
nians and the NCSM3 (triangles) NN starting Hamiltonian, all 
for N max =6. 
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FIG. 6. (Color online) Calculated exclusive cross sections 
(upper panel) and FWHM widths of the momentum distribu- 
tions (lower panel) for the ground and 7-decaying final states 
of the 10 B residues, following np-removal at 2100 MeV per 
nucleon. The TNAs used were obtained using the WBP in- 
teraction (solid lines and solid squares), the pjt interaction 
(solid lines and open squares) shell-model interactions and the 
NCSMl (dashed lines blue circles) and NCSM2 (dashed lines and 
open circles) NN+3N Hamiltonians for N ma x=&- 



ble IIV1 the simplest case being the 3 + ground state, 
where only a single p-shell configuration contributes. The 
[OP3/2] 2 TNA vary, depending on the different interac- 
tions, and the magnitude of these p-shcll amplitudes is a 
key factor in the cross section changes observed. 

For the two T = 0, 1 + states, a mix of p-shell config- 
urations now contribute, with the overall magnitude of 
the TNA and their relative strengths and phases chang- 
ing. The relative strengths of the three configurations 
are broadly consistent across the NCSM interactions, but 
are different from the truncated-basis WBP interaction. 
In particular the [P3/2] 2 TNA arc different, with some 
apparent shift of strength from the first to the second 1 + 
state, when compared to the NCSM interactions. In these 
cases there is interference between the different configura- 
tions that make it less transparent what one expects from 
the different interactions. Changes in the p-shell config- 
urations will account for a part of the changes in cross 
sections shown in Table IIIII and the predicted momen- 
tum distribution widths for the second 1 + state, offering 



a means to discern between the different interactions. 

Further enhancement of the cross section may arise 
from the new couplings to higher major shells. Cou- 
pling to major shells of the same parity (odd N) will lead 
to generic changes to the overall size of the two-nucleon 
overlap functions and TNA. The TNA due to mixing with 
major shells of opposite parity (even N) leads to new in- 
terference affects that can enhance two-nucleon spatial 
correlations (see e.g., |16|). For the first 1 + this can be 
seen in Fig. [7]for calculations based on the NCSMl Hamil- 
tonian. The left panel is calculated when retaining only 
the Qhaj p-shcll TNA components from the N max = 6 
NCSM calculation. The right panel includes the full set 
of NCSM TNA for all major shells. The enhanced spa- 
tial correlations presented to the target nucleus from the 
inclusion of single-particle configurations with opposite 
parity in the two-nucleon overlap function are evident. 
The cross sections from these truncated and full TNA 
sets are 6.09 and 8.52 mb, respectively. Both exceed 
those of the p-shell shell-model calculations, this being 
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TABLE III. Like and unlike two-nucleon removal cross sec- 
tions (in mb) for a C projectile incident on a carbon target 
at 2100 MeV per nucleon. The excitation energies, Ef, of 
each final state are shown in Fig. [T] The TNAs used were 
calculated using (a) the Ofku p-shell-model and wbp inter- 
action, (b) the two NN+3N interactions NCSMl and NCSM2 
(see text), and (c) the NCSM3 Hamiltonian in which the chi- 
ral 3N interaction is turned off. Calculations (b) and (c) use 
the NCSM with N max = 6. The sums show the accumulated 
cross sections that lead to the ground state and the 7-decaying 
bound excited states of the mass A — 10 projectile residues. 
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state decays by a emission with a 16% a-branch. 



TABLE IV. p-shell two-nucleon amplitudes for the wbp, 
NCSMl, NCSM2 and NCSM3 interactions. 



Interaction 


(J" 


T) 


[pi/2] 2 


[Pl/2][P3/2] 


ba/2] 2 


WBP 


(3f 


0) 






1.976 




(If 


0) 


-0.011 


0.979 


0.699 




(If 


0) 


0.363 


0.229 


-1.134 


NCSMl 


(3f 


0) 






1.913 




(If 


0) 


-0.220 


1.034 


1.197 
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0) 


0.611 


0.376 


-0.835 


NCSM2 


(3f 
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2.213 




(If 
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-0.255 


0.863 


1.307 




(If 
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0.470 


0.500 


-0.814 


NCSM3 


(3f 
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1.644 




(If 


0) 


-0.224 


1.137 


1.205 




(If 


0) 


0.740 


0.332 


-0.719 



3.58 mb for the WBP interaction TNA. 

For both NCSMl and NCSM2 the enhancement of the 
np removal cross sections relative to the wbp calcula- 
tions is larger than that for the like-nucleon, nn, removal 
cross section. However, the difference in this ratio is rel- 
atively small. Despite the larger np removal cross section 
obtained, some underestimation of the experimental np 
channel cross sections still remains. The available data 
again suggests that there are remaining deficiencies in 
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FIG. 7. (Color online) Projected two-nucleon density 
Vf(si, S2) from the two-nucleon overlap for the first T = 0, 
1 + state, calculated using the NCSM and the NCSMl Hamil- 
tonian. The contours show the position probability in the 
impact parameter plane of nucleon 2, S2, when nucleon 1, si, 
is at the position indicated by the black point. The left panel 
is calculated when retaining only the lowest p-shell TNA from 
the Nmax = 6 calculation. The right panel includes the full 
set of NCSM TNA for N max = 6. The overlap shows the 
enhanced spatial correlations arising from inclusion of single- 
particle configurations of opposite parity in the two-nucleon 
overlap function (see also Fig. 4 of Ref. The cross 

sections from these two TNA sets are 6.09 and 8.52 mb, re- 
spectively. 



the T = parts of the two-nucleon overlap functions. 
The yields to specific final states, namely the first 3 + , 1 + 
and + states, suggest a significant sensitivity to the in- 
teractions used, requiring input from final-state-exclusivc 
measurements . 



V. SUMMARY 

We have considered the impact of microscopic, NCSM 
wave function overlaps on the theoretical cross sections 
for two-nucleon removal reactions from fast 12 C projec- 
tiles. Data were available for reactions on a carbon target 
at beam energies of 250, 1050 and 2100 MeV per nucleon. 
As found in a previous analysis @ , the np removal cross 
sections are underestimated by the theoretical model cal- 
culations, but do show an enhancement relative to the use 
of truncated-basis p-shell-modcl calculations. The cross 
sections to both T = and T = 1 states are enhanced, 
and the use of large-basis NCSM amplitudes does not fully 
resolve the relative discrepancy between measured np- 
and nn-removal cross sections. 

Further measurements, of final-state-exclusive cross 
sections and residue momentum distributions, would al- 
low a much more detailed scrutiny and confrontation of 
the detailed reaction observables predictions, including 
the identification of any indirect reaction components 
arising from two-step paths to the final-states. The cal- 
culated np-removal cross sections to the T=0, 10 B final 
states were shown to have sensitivity to the different vari- 
ants of the chiral interactions used; for example, the ra- 
tio of the calculated cross sections to the 10 B ground 
3 + , T = and first 1 + . T = excited states. To a lesser 
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degree the first + , T = 1 state cross section and the 
branching between the first T = 1, + and 2 + states was 
found to depend on the effective interaction. In this case, 
data for the nn and pp removal channels would provide 
useful verification. The momentum distribution of the 
second 1 + , T = state also shows a particular sensitivity 
to the interaction, providing a further useful probe. 

The overall conclusion from the present analysis is 
that the existing rcsiduc-final-state-inclusive data sug- 
gest that the T=0, np-spatial correlations present in the 
wave functions used are still insufficient. We have shown 
that new exclusive measurements would offer a means 



to interrogate these shell-model inputs, in particular for 
the np-channcl, T = wave functions, and the direct 
reaction mechanism predictions in considerable detail. 
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